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Among vertebrate gastrointestinal microbiome studies, complete representation of taxa is limited, 
particularly among reptiles. Here, we provide evidence for previously unrecognized host-microbiome 
associations along the gastrointestinal tract from the American alligator, a crown archosaur with shared 
ancestry to extinct taxa, including dinosaurs. Microbiome compositional variations reveal that the digestive 
system consists of multiple, longitudinally heterogeneous microbiomes that strongly correlate to specific 
gastrointestinal tract organs, regardless of rearing histories or feeding status. A core alligator gut 
microbiome comprised of Fusobacteria, but depleted in Bacteroidetes and Proteobacteria common to 
mammalians, is compositionally unique from other vertebrate gut microbiomes, including other reptiles, 
fish, and herbivorous and carnivorous mammals. As such, modern alligator gut microbiomes advance our 
understanding of archosaur gut microbiome evolution, particularly if conserved host ecology has retained 
archosaur-specific symbioses over geologic time. 

Vertebrate gastrointestinal (GI) tracts are distinct microenvironments formed by conserved ecological 
interactions and ancient evolutionary processes between a host and resident microbial flora. Interest in 
unraveling the ecologies and evolutionary histories of host-microbe symbiotic associations has surged 
following the advent of next-generation sequencing and '-omics' approaches 14 . However, reconstruction of 
vertebrate-microbe evolutionary histories has been limited because some host groups are underrepresented 
among microbiome studies. Specifically, despite mammals representing only —10% of known vertebrate taxa 5 , 
the majority of GI tract microbiome studies to date focus on mammalian, predominately human, health and 
disease (e.g. obesity) 6 or microbiome acquisition 7 . In contrast, reptilian taxa (e.g. crocodylians, tuatara, turtles, 
squamates) 810 represent —17% of known, extant vertebrate species 5 . Except for evaluation of fecal coliform 
bacteria 11 , which typically do not reflect the actual GI tract microbial diversity 12 , and from culture-based studies 
focusing on Salmonella spp. 13 , almost nothing is known about modern crocodylian GI tract microbiome diversity 
and ecology, despite its importance 14 to advance respiratory and cardiovascular system development studies 15 . 

Here, we expand the current knowledge of GI tract microbiomes from the semi-aquatic, predominantly 
freshwater American alligator (Alligator mississippiensis), indigenous to the southeastern United States. From 
wild and farm-raised alligators from fasting and feeding seasons, we identify distinct bacterial communities 
partitioned as a function of organ type along the length of the GI tract and the presence of a shared or "core" 
bacterial community. Because the modern alligator has been used as an analogue to understand and infer 
physiological 16 and biomechanical 17 traits of extinct aquatic and terrestrial forms, including carnivorous theropod 
dinosaurs 17 , as well as to characterize extinct archosaur thermoregulation 18 , feeding 19 , respiration 16 , and mobil- 
ity 17 , we propose that the alligator microbiome also provides insight into the evolutionary history of predator 
digestive systems 14,20,21 and archosaur microbiome symbioses 22 based on conserved ecological niche and extant 
symbiotic associations 2,23,24 . 

Results 

Wild A. mississippiensis generally fast from October to March 25 . We hypothesized that protracted fasting would 
yield the composition of the indigenous symbiotic gut microflora, based on findings from the Burmese python 
that captured postprandial remodeling in microbiome composition 8 . During fasting 25 , the wild alligator stomach 
contents (<15 mL) consisted of viscous, semi-opaque yellow gastric juices with pH 2.1 to 3.0 (average 2.6; n = 5) 
and that were occasionally mixed with green, fluid-like material reminiscent of vegetation. Alligators routinely 
ingest vegetation as a consequence of capturing prey along shorelines, and prey for wild alligators includes small 
mammals, birds, crustaceans, fish, amphibians, and other reptiles 26 . During feeding, wild alligator stomach 
contents consisted of residual plant material, cutaneous fragments of partially digested crawfish, and whole 
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crawfish, with a stomach pH 1.6 to 3.7 (average 2.8; n = 10). From 
animal "N", there were 17 whole, 2.5 cm long crawfish recovered, 
and from animal "D", nematodes, common in wild alligators, were 
recovered. In contrast, the farm-raised animals were fed a commer- 
cially available, dry, pelleted alligator ration (manufactured by 
Cargill) that consisted primarily of crude protein (45-56%), fat (9- 
12%), and fiber (3-4%) derived from animal protein products (dried 
blood meal, meat, and bone), wheat, soy, corn, and animal fat. 
Stomach contents from farm-raised individuals were viscous, opa- 
que, brown to yellow fluids, and stomach pH ranged from 1.7 to 6.6 
(average 3.3; n = 5) during fasting and from 3.0 to 5.9 (average 4.4; 
n = 5) while feeding. 

General patterns of bacterial phylum-level and proteobacterial 
class-level relative abundances for each alligator microbiome corre- 
lated to seasonal fasting or feeding based on nonparametric multi- 
variate analysis of variation (NPMANOVA) (Fig. 1; Supplementary 
Tables 1 and 2). Wild fasting (winter) microbiome compositions 
shifted from communities dominated by Proteobacteria to Firmi- 
cutes- and Fusobacteria-dominated communities after renewed 
feeding in the spring (NPMANOVA P = 0.001, F = 8.17). 
Gammaproteobacteria comprised —30% of all wild (fasting) tissue 



communities, whereas Firmicutes and Fusobacteria comprised 67% 
and 25%, respectively, for wild spring-collected tissues. Firmicutes 
represented 17.6% of the winter wild microbiome amplicons, 
but 66.9% upon renewed feeding and active fat deposition 27 
(NPMANOVA P < 0.001, F = 16.75). In contrast, Firmicutes com- 
prised a significant portion of both winter farm-raised (35.5% of 
recovered amplicons) and spring farm-raised (48.5% of recovered 
amplicons) gut communities (NPMANOVA P = 0.002, F = 3.66). 
Animal rearing history (i.e. wild or farmed), sex, and overall length 
(as a proxy for age and size) did not significantly explain the observed 
bacterial variability (NPMANOVA P > 0.05). Similarly, non-metric 
multidimensional scaling (NMDS) plots showed that wild and 
farm-raised bacterial communities responded similarly to scaling 
dimensions in NMDS space (Supplementary Fig. 1). Among all the 
alligators, seasonal feeding/fasting status explained the significant 
differences in bacterial compositions (NPMANOVA P < 0.001, 
F = 6.09). 

But, for each alligator, significant changes in bacterial community 
compositions were attributed to specific organ or tissue type (i.e. in 
reference to changes in epithelial tissue along the alimentary canal) 
(NPMANOVA P < 0.001, F = 2.42; Fig. 1). The mouths had the 
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Figure 1 | Phylum-level, gastrointestinal bacterial community representation (>2% of retrieved amplicons) from wild and farm-raised, winter 
and spring, A. mississippiensis individuals and schematic sample locations. Data are divided based on rearing history (wild vs. farmed), and 
subdivided based on season (winter and spring, blue and pink backgrounds, respectively) (Supplementary Table 1). Significant differences in bacterial 
composition between downstream tissues or organs (paired Student's f-test) are denoted by asterisks: P = 0.05 — 0.01 (*), 0.01 — 0.001 (**), and P < 
0.001 (***) (full f-tests in Supplementary Table 7). 



SCIENTIFIC REPORTS | 3 : 2877 | DOI: 10.1038/srep02877 



2 



richest tx-diversity based on the number of operational taxonomic 
units (OTUs). As alligators frequently open their mouths for ther- 
moregulation 28 , rich a-diversity potentially reflects frequent inocu- 
lation with transient environmental bacteria. The upper GI tract 
communities (e.g. stomach tissue, fluids, duodenum) had the lowest 
richness and diversity (Supplementary Fig. 2), likely as a con- 
sequence of low pH. Intermediate richness and a-diversity values 
were obtained for the lower GI tract tissues (colon and feces; 
Supplementary Fig. 1). From the oral and upper GI tract only (i.e. 
stomach, duodenum), 94% of the amplicons were affiliated with 
Proteobacteria (Betaproteobacteria and Gammaproteobacteria), 
Bacteroidetes, and Firmicutes, whereas 73% of the amplicons from 
the lower GI tract (i.e. ileum, colon, and feces when present) were 
affiliated with Fusobacteria, Firmicutes, and Bacteroidetes. In gen- 
eral, Fusobacteria dominated the lower GI tract, particularly feces 
(NPMANOVA P = <0.001, F = 23.5), but were a minor group (4%) 
when present in the rest of the GI system. Unique community 
ordination according to tissue type in NMDS space (Fig. 2) showed 
that communities retrieved from the colon and ileum occupied the 
largest NMDS ordination space compared to other tissue communit- 
ies that did not overlap (i.e. feces and tongue). This was supported by 
(3-diversity analyses of shared bacterial OTUs that also revealed a 
potentially high level of endemism to each organ. The greatest sim- 
ilarity levels in community composition were only among physio- 
logically adjacent organs or from samples derived from the same 
organ (i.e. stomach and gastric juices; Supplementary Table 3). 

We used average-neighbor cluster analyses of amplicons recov- 
ered from each organ or tissue type to assess the potential for a core 
microbiome. Defining a host-specific or environment-specific core 
microbial community can be used to develop a screening tool that 
could identify host health, immunity, and disease 29,30 and to deduce 
symbiotic evolutionary relationships through time 31 . The number of 
unique OTUs were not significantly different between feeding and 
fasting in wild (two-tailed, paired r-test P = 0.17) and farmed (P = 
0.25) animals. More unique OTUs were recovered from oral samples; 
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Figure 2 | Euclidean-based NMDS plot of bacterial diversity as a function 
of tissue type from all individuals, regardless of rearing history or feeding 
status. Gastrointestinal samples proximal to each other had more similar 
bacterial communities. The plot was constructed using normalized log 
abundances, and the parameters season (spring or winter), animal length, 
rearing history (farmed or wild), and sex (male or female). Lines denote 
vector overlays on NMDS ordination, indicating the directionality and 
strength of change (line length) for a specific parameter. 



however, for individuals where tongue tissues were unavailable and 
esophageal tissues were utilized, the greatest number of unique 
OTUs were from the lower GI (colon or feces; Supplementary 
Table 4). The number of unique OTUs specific to wild and farm- 
raised organ or tissue type ranged from 2.7% up to 86.6% of all 
recovered OTUs from that sample type (average of 26.2% unique 
OTUs per organ or tissue; Supplementary Table 3). The largest num- 
ber of shared OTUs based on average-neighbor cluster analysis was 
observed in individuals from the same season (e.g. winter vs. spring; 
Supplementary Table 4). A maximum of 13 OTUs were shared 
among all wild individuals (Fig. 3; Supplementary Table 5). 
Although it may be too stringent given potential error incurred with 
pyrosequencing 32 , a 96% sequence similarity cutoff was used to 
denote species- to genus-level associations 33 among the OTUs. One 
OTU from the GI tracts of all wild and farmed individuals was 
assigned to the Enterobacteriales. The number of shared OTUs 
increased by lowering the sequence similarity cut-offs to 95% and 
91%, which represented higher taxonomic associations 32,33 (Table 1). 
Resulting OTUs identified at this level included representatives of 
Bacteroidetes (Bacteroidales), Firmicutes (Clostridia, Clostridiaceae, 
Clostridiales), Fusobacteria (Fusobacteriaceae), and Gammaproteo- 
bacteria (Enterobacteriales). 

Discussion 

Microbiome acquisition occurs over evolutionary time and reflects 
complex feedbacks between the host 34 , environment 35 , diet 36 , im- 
mune response, and microbe-microbe community interactions 37 that 
change as new lineages invade new ecological niches. However, rep- 
tilian taxa are currently underrepresented among gut microbiome 
studies, thereby resulting in limited potential to reconstruct the 
evolutionary history of vertebrate-microbe symbioses. Based on pre- 
vious studies that emphasize host feeding status controls microbiome 
composition 8 , and because alligators have distinctive physiological 
and biochemical differences compared to other animal hosts, par- 
ticularly with respect to seasonal fasting 38 , we expected that alligator 
GI tract microbiomes would differ as a function of rearing history. 
Our results reveal longitudinally heterogeneous microbiomes along 
the alligator GI tract, with compositional differences significantly 
correlating to organ or tissue and also to the farm-raised versus wild 
alligator diets. 

The dissimilarity between our alligator microbiome compositional 
findings and the conclusions drawn in earlier GI tract studies is 
striking. This is likely because the previous work, including from 
mammalian hosts, relies predominately on colonic biopsies 12,29 31 , 
feces 9,23,31,39 , gastric fluids 40 , cecal contents 3 , or intestines 1,29,35 to 
define host microbiomes. However, these limited samples may not 
provide the insight intended 12 . From our analyses, distinct microbial 
communities were detected between the alligator mucosa and lume- 
nal contents, as has also been observed in mice 41 . Also, the fecal 
microbial community composition did not represent the composite 
alligator GI tract microbiome. Representatives from the Fusobacteria 
phylum were most abundant from farm-raised and wild feces, but 
were minor members of GI tract communities (Supplementary Table 
1). In contrast, other undetected phyla in the feces (e.g. Actinobac- 
teria) had significant representation in other regions of the GI tract. If 
alligator feces alone were used to represent the gut microbiome, or at 
least the lower GI tract, then Bacteroidetes would be considered the 
most prevalent microbial group in the alligator microbiome. But, 
four OTUs total were affiliated with Bacteroidetes, which comprised 
only 10% of the composite microbiome from all of the digestive 
organs (Table 1). The relative abundances of Bacteroidetes were 
considerably lower than what has been reported in other vertebrate 
microbiome studies 40,42 , which again may be due to the use of fecal, 
cecal, or intestinal samples alone. Additionally, it is possible that 
increasing the number of individual hosts in a study could alter 
interpretations of inferred composite microbiomes. However, at least 
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Figure 3 | Schematic average neighbor cluster analyses for shared 
bacterial OTUs retrieved from wild and farm-raised alligators. OTU 

taxonomy was based on BLAST results (Supplementary Table 5), and 
the value shown in brackets refers to the number of OTUs at phylum- 
or class-level. Cluster analysis at 96% sequence similarity (solid black 
lines) resulted in 63 shared OTUs, one of which was recovered from all 
wild and farmed individuals (Enterobacteriales). Reducing the similarity 
threshold to 95% (dashed lines), and 91% (grey lines), reduced the 
number of OTUs to 59 and 41, respectively. At 95%, two shared OTUs 
among all animals formed. 



from our results, the observed pattern of microbial community vari- 
ation along the length of the alimentary canal should remain because 
of the underlying function of each organ for digestion and the sig- 
nificant correlation between community composition and organ or 
tissue type. Future studies should utilize representative samples of 
the entire gut if the GI tract is considered to be a heterogeneous 
ecosystem with organ-specific microbiomes. 

Significant differences in the alligator GI tract microbiome com- 
munity composition according to organ or tissue type may be due to 
the host's metabolism and physiology (e.g. carnivore, herbivore, or 
omnivore) and metabolic capabilities of members of the gut micro- 
biome. But, differences between farm-raised and wild alligator gut 
bacterial community compositions could also be attributed to host 
diet. Firmicutes dominated the farm-raised alligator GI tract micro- 
biomes throughout the year, comparable to microbiomes of other 
captive-bred and raised reptiles 8,9 . However, Firmicutes only domi- 
nated wild alligator microbiomes during seasonal feeding. The farm- 
raised alligators had significantly greater fat deposits compared to 
their wild counterparts. Obesity in farm-raised alligators has been 
observed previously 26 due to a diet high in carbohydrates and satu- 
rated fats 38 to increase animal size for commercial meat and skin 
production. Obesity and excess fatty tissue have been correlated to 
shifts in microbiome composition in humans 6,42 . In other mam- 
malian taxa, including mice 42 , Firmicutes are also prevalent members 
of microbiome communities from obese individuals. The intriguing 
association between increased fat deposition and the predominance 
of Firmicutes in gut microbiomes from feeding alligators, as well as 
from obese vertebrates, suggests that diet has the potential to modify 
gut microbiome physiology when the ecological niche occupied by 
the wild counterpart changes. 

From the alligator GI tract microbiomes, a core community of 
Fusobacteria, Bacteroidetes, Firmicutes, and Proteobacteria is strik- 
ingly different from the core microbiomes inferred for other verte- 
brates, including mammals 23,43 , birds 3-4,39 , and even other reptiles 8-10 
(Fig. 4, Supplementary Tables 7, 8). As all other vertebrates sampled 
to date possess representatives of Bacteroidetes, Firmicutes, and 
Proteobacteria, regardless of host taxonomy, the presence of 
Fusobacteria in all alligators sampled makes the alligator gut micro- 
biome unique. Fusobacteria represented —12% of the OTUs in the 
core alligator microbiome, and dominated the lower GI tract tissue 
bacterial communities and fecal material (Table 1). In contrast, 
Fusobacteria have been retrieved in low abundances from composite 
gut microbiomes from human (<10% of all oral sequences) 43 , fish 
(<13% of sequences) 35 , and chicken (<3.0% of sequences) 4 hosts, 
with Fusobacteria recovered from human oral cavities playing a crit- 
ical role in initial biofilm development 44 . In these other hosts, 
Firmicutes and Bacteroidetes represent the overwhelming majority 
of recovered microbiome OTUs. 

Our results raise a fundamental question: why are Fusobacteria 
prevalent in the American alligator GI tract? Fusobacteria are a 
poorly-studied phylum comprised of approximately 32 described 
species, with an overall uncertain phylogenetic history. Some 
researchers place the phylum at a basal position 45 based on both 
rRNA and core protein analyses with the hyperthermophilic Aqui- 
ficae and Thermotogae phyla, but others consider Fusobacteria 
to have closer phylogenetic affiliation with Bacteroidetes and 
Firmicutes 44 . Placement as a basal phylum has led to an inferred 
divergence —3.5 billion years ago 45 . However, because of high hori- 
zontal gene transfer within Firmicutes, and because Fusobacteria are 
frequently recovered with Firmicutes, a divergence no earlier than 
400 Myr ago is also proposed, which coincides with the origin of the 
vertebrate digestive tract 44 . The almost complete lack of envir- 
onmental studies that describe the presence, diversity, and functional 
role of Fusobacteria 45 in microbiomes confounds the uncertain 
phylogenetic position. Although Fusobacteria are associated with oral 
diseases in humans 44 , their prevalence in alligators can be interpreted 
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Table 1 | Taxonomic assignments of shared OTUs 
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The numbers of OTUs refer to the total number of OTUs recovered in wild or farm-raised individuals and assigned to the phyla or class, where assignment was possible according to BLAST (Supplementary 
Table 5). Values in parentheses refer to the percentage of the shared community represented by each phylum. 



as specialized, perhaps critical to host health or nutrient acquisition, 
because alligators possess innate immune compatibility not present 
in mammalians. Fusobacteria are also known to play a role in biofilm 
development 44 , so their presence in the GI tract could affect lumen 
biofilm development. This would explain the dominance of Fuso- 
bacteria in alligator feces. As other vertebrate microbiome studies 
rarely sample epithelial tissue (although see 46,47 for exceptions), it is 
possible that Fusobacteria have been undersampled in microbiome 
studies until now. Nevertheless, based on what data are available, 
we speculate that Fusobacteria in the lower GI tract of alligators 
may occupy a functional role in digestive organ development and 
nutrient acquisition that precedes a similar ecological niche that 
is now occupied by Firmicutes and Bacteroidetes in mammals based 
on a combination of strong selective pressures driven by host gen- 
etics, a conserved diet, niche occupation, and unique biochemical 
adaptations 48 . 

Microbiomes have the potential to provide a window into an 
ancestral condition that has been retained over geologic time. 
Others postulate that vertebrate gut microbiome evolution may cor- 
relate to tetrapod evolution 2,23 . But, the diversity of bacterial com- 
munities residing in vertebrate guts is far from explored, as evidenced 
from our study that reveals a distinct alligator gut microbiome. An 
analysis of current microbiome data, albeit based on feces that likely 
represent biased microbiome compositional information (Fig. 4), 
still shows that microbiome compositions do not reflect host phylo- 
genetic affiliations for all Vertebrata. Placement of sister taxa based 
on microbiome composition alone may lead to non-biologically 
meaningful information (i.e. placement of gorillas as sister taxa to 
zebrafish rather than to other hominids included in the dataset). 
Consequently, what is known from mammals is likely not represent- 
ative of tetrapods as a whole, although successful attempts have been 
made to reconstruct the evolutionary history of mammalian digestive 
systems 2,23,24 . 

With a similar approach, the alligator microbiomes can be used 
conservatively to reconstruct the ecology and evolutionary history of 
crocodylian, and perhaps even archosaur, digestive systems 20 . 
Pseudosuchia and Avemetatarsalia, lineages within Archosauria, 
diverged in the Early Triassic (approx. 250 Myr ago), giving rise 
to modern crocodylians and avians, respectively 49 . The order 
Crocodylia is comprised of 23 extant species of caiman, alligators, 



crocodiles, and gharials within the Alligatoridea, Crocodyloidea, and 
Gavialoidea families 50 . Recent phylogenetic reconstructions from 
extant taxa and the robust crocodylian fossil record, which extends 
from the Late Cretaceous (70-84 Myr ago), highlight that modern 
crocodylians have had extensive, dynamic evolutionary histories that 
inform us about the ecology, physiology, and biochemistry of the 
ancient lineages 18,50 . For instance, the biochemically unique croco- 
dylian blood serum with its antibacterial properties was likely an 
ancestral trait that originated in response to aggressive behavior 
associated with frequent injuries and occupation of potentially 
pathogen-rich wetlands 48 . Fossil evidence indicates that numerous 
crocodylian taxa had similar trophic statuses and occupied compar- 
able ecological niches to the modern 51 , generally functioning as top 
predators and scavengers within a narrow aquatic to semi-aquatic 
niche (e.g. wetlands and floodplains). We do not know when modern 
archosaur symbioses became established in the geologic past, and 
what the ancestral condition may have been. However, it is possible 
that ancestral crocodylian and extinct archosaur GI tract micro- 
biomes would have been more similar to modern alligators than to 
other vertebrates studied to date, assuming comparable trophic sta- 
tus, diet, and niche occupation through time. Future research should 
include a robust evaluation of the ecological status (i.e. aquatic, semi- 
aquatic or terrestrial), feeding habits, and phylogenetic position 
of archosaurs 49,52 to test this hypothesis. Nevertheless, the implica- 
tions from the modern alligator gut microbiome comprised of 
Bacteroidetes, Firmicutes, Proteobacteria, and Fusobacteria, the 
dominance of which is unique among vertebrates, advance our 
understanding of archosaur gut microbiome evolution and host- 
microbe symbioses. 

Methods 

Sample and dataset acquisition. Organs and tissues sampled along the length of the 
gastrointestinal tract (Fig. 1) from eight salvaged A. mississippiensis juveniles (114— 
180 cm, snout to tail lengths) were acquired within hours of death from individuals 
sacrificed for other studies in December 2009, May/ April, August, and December 
2010 at the Rockefeller Wildlife Refuge, Louisiana, USA (under the Louisiana Natural 
Heritage Program, permit LNHP-10-009). In the winter, two farm-raised but fasted 
(140 and 160 cm, snout to tail length; fasted 6 weeks) and two wild fasting (114 and 
180 cm) alligators were used. In the spring, two fed farm-raised (130 and 131 cm) and 
two presumed feeding (which was later confirmed from the presence of prey items in 
the stomachs) wild alligators (132 and 140 cm) were used. As the alligators were 
salvaged, not all tissues were available for sampling every time (e.g. esophageal tissue 
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Figure 4 | Euclidean -based cluster analysis of phylum-level bacterial 
diversity of representative vertebrate taxa. Colors correspond to major 
divisions of host phylogenetic affiliation (mammals, reptiles, birds, and 
fish) and letters reflect diet. Data used for alligator are based on average 
feces bacterial composition (Supplementary Table 8). Feces were selected 
as the majority of microbiome studies utilize fecal material. Phylum-level 
bacterial assignments were selected as the published data generally did not 
include lower taxonomic assignments. 



was collected from a wild and a farmed animal during the spring because tongue 
tissue was unavailable due to skull removal by a different researcher). Within hours of 
death, several grams of epithelial tissue were aseptically excised or scraped from along 
the digestive tract with sterile scalpel blades and stored in 100% molecular grade 
ethanol at ~20°C. Tissue scrapings ensured complete evaluation of the biofilms lining 
the GI tract, as well as any microbes adhering to or within the mucosa. Gastric fluid 
was collected by making a small (6-10 cm length) incision into the stomach with 
sterile scalpel blades to provide direct access to the contents. Liquids were collected 
using sterile pipettes and solid material was collected using sterilized forceps. 
Stomach contents were then transferred to sterile 50 mL Falcon tubes. Following the 
collection of gastric contents, tissue scrapings were collected. Fecal samples were 
collected directly from the colon aseptically with sterile scalpel blades by making a 
small (6-10 cm) incision into the colon. As with gastric tissue samples, colon 
scrapings were collected following collection and removal of fecal material. In 
triplicate, tissue (~1 cm 3 ), gastric fluids (200 ul), and fecal material ( — 1 cm 3 when 



present) were homogenized in 150 to 200 ul lx TE buffer using a hand-held Kontes 
pellet pestle, and total nucleic acids were extracted using the Qiagen DNeasy kit, 
following manufacturer instructions. 

Purified DNA from each gastrointestinal sample (n — 48) was sent to the Research 
and Testing Laboratories in Lubbock, TX (USA), for 454 tag-encoded FLX titanium 
amplicon pyro sequencing of the V1-V3 region of bacterial 16S rRNA genes from 
nucleic acids extracted from gastric fluids, fecal material, and epithelial lining from 
along the digestive tract. Samples from two individuals (SKA_09D_W and 
SKA_09D_F) were sequenced in January 2010 and the remaining samples from six 
additional individuals were sequenced in January/February 2011. After purification 
of DNA using previously described methods 53 , amplicon pyro sequencing was done 
using the forward primer 28f (5 ' -GAGTTTGATCNTGGCTCAG-3 ' ) and the reverse 
primer 519r (5' -GTNTTACNGCGGCKGCTC-3 ' ) 53 . 

Dataset processing and statistical analysis. Amplicon trimming, removal of 
primers, barcodes, and low quality reads, screening for chimera, alignments, 
taxonomic assignments, assessment of a- and ^-diversity, and statistical analyses 
were consistent with previous studies 3,30,54 ' 55 . Additional screening, OTU definitions 
at 96% sequence identity, average- neighbor cluster analysis, rarefaction curves, and 
statistical calculations were run with mothur, version 1.28.0 54 . The datasets were used 
to calculate rarefaction curves (Supplementary Fig. 3), revealing high to moderate 
sample coverage, as well as bacterial diversity indices (Supplementary Fig. 2) and 
richness values (Supplementary Fig. 2). The Ribosomal Database Project (RDP) was 
utilized for sequence trimming, alignment, and taxonomic assignment 3,30 . Processing 
included screening to trim reads less than 150 bp in length, removing primers, 
barcodes, and low quality reads (Q < 20) using RDP. Following alignment in RDP, 
potentially chimeric reads were identified using the uchime command in mothur 54 . 
Reads identified as potentially chimeric were visually examined, compared to known 
reads using BLAST, and removed. Additional screening following chimera removal 
using mothur included removing sequences with homopolymeric regions (>8 bp), 
removing ambiguous reads, and ensuring uniform start positions for all reads and a 
minimum length of 150 bp 55 . After filtering, reads were uploaded to RDP for 
taxonomic assignment. Subsequent average -neighbor clustering, rarefaction curve 
generation, and statistical analyses were performed using mothur. A core 
microbiome was obtained by clustering reads using variable OTU cut-offs (96%, 95%, 
91% sequence similarity). Raw pyrosequence files were deposited in the GenBank 
Short Read Archive (accession numbers SRA023831 and SRA062824). Student's t- 
tests and NPMANOVA were performed in PAST 56 to test the significance of bacterial 
community composition, number of recovered OTUs, and percentages of OTUs (a = 
0.05). NPMANOVA based on Euclidean distances of log-transformed data was 
utilized to evaluate the environmental variable(s) that accounted for the observed 
variations in bacterial community composition. 
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